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Abstract 

The responses of long-term growth of plants under elevated C0 2 have been studied extensively. Comparatively, the 
responses of plants to subambient C0 2 concentrations have not been well studied. This study aims to investigate the 
responses of the model C 3 plant, Arabidopsis thaliana, to low C0 2 at the molecular level. Results showed that low 
C0 2 dramatically decreased biomass productivity, together with delayed flowering and increased stomatal density. 
Furthermore, alteration of thylakoid stacking in both bundle sheath and mesophyll cells, upregulation of PEPC and 
PEPC-K together with altered expression of a number of regulators known involved in photosynthesis development 
were observed. These responses to low C0 2 are discussed with regard to the fitness of C 3 plants under low C0 2 . This 
work also briefly discusses the relevance of the data to C 4 photosynthesis evolution. 
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Introduction 

The response of plants grown in lower C0 2 concentrations 
has been much less studied than responses to elevated C0 2 
concentrations (Long et al. , 2004, 2006; Ainsworth and Long, 
2005; Gerhart and Ward, 2010). Among these limited stud- 
ies, some have demonstrated that a large genetic variation in 
response to low C0 2 exists among Arabidopsis accessions. For 
example, Sharma et al. (1979) screened 33 Arabidopsis acces- 
sions for survival time under limiting C0 2 when grown side 
by side with C 4 plants (Zea mays L.) in an air-tight chamber 
where C0 2 concentration was reduced to below the compen- 
sation point of C 3 plants and found a 1-2-week difference in 
the survival time in different accessions and also found sub- 
stantial genetic segregation among F 2 parents, with extreme 
differences in survival time near the C0 2 compensation point. 

Arabidopsis genotypes from different elevations show sig- 
nificant variation in the response of seed number when grown 
at low CQ 2 (20 Pa) (Ward and Strain, 1997). Ward et al. (2000) 
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performed an artificial selection experiment using Arabidopsis 
for high seed number over five generations at low C0 2 (20 Pa, 
or 200 ppm); the selected populations produced 25% more 
seeds and 35% more biomass on average than control popu- 
lations which were randomly selected at the fifth generation 
when grown at low C0 2 . In addition, Ward and Kelly (2004) 
also observed a high level of genetic variation in survival, 
reproductive output, and total seed production among the 
Arabidopsis genotypes when grown at low C0 2 (200 ppm). 
All these studies suggest that Arabidopsis has adaptive phe- 
notypic plasticity in response to low C0 2 . 

In a carbon starvation experiment, 5-week-old Arabidopsis 
rosettes treated with ambient (350 ppm) C0 2 or compensa- 
tion point (<50 ppm) C0 2 were collected in the light for 4h 
to investigate responses to changing endogenous sugar con- 
centrations in rosettes at the gene expression level using the 
GeneChip Arabidopsis ATH1 genome array (Biasing et al, 
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2005). However, these studies have not addressed the mecha- 
nism of long-term responses of plants to low C0 2 . 

This study conducted a survey of responses of C 3 plants 
to long-term low C0 2 treatments at the molecular level. 
Arabidopsis was chosen as the model system because its 
genome has been fully sequenced and is still the best annotated 
plant genome to date (The Arabidopsis Genome Initiative, 
2000); the well-annotated Arabidopsis genome facilitates 
analysis of global gene expression using RNA-Seq technol- 
ogy. This study sequenced the transcriptome of 6-week old 
Arabidopsis seedlings grown under ambient C0 2 (380 ppm) 
or low C0 2 (100 ppm). The results are discussed with particu- 
lar reference to the significance of the altered gene expression 
to the fitness of C 3 plants under low C0 2 . The relevance of 
low C0 2 to C 4 evolution is also briefly discussed. 

Materials and methods 

Plant growth and harvest 

Arabidopsis thaliana Columbia-0 (Col-0) seeds were imbibed in 
0.1% (w/v) agar solution and incubated at 4 °C for 2 d to break dor- 
mancy. Imbibed seeds were germinated and grown in Pindstrup soil 
in a Percival incubator (NC-350HC-LC, Nihonika, Japan) in which 
C0 2 gas can be accurately and stably controlled in the range of 
100-3000 ppm. C0 2 concentrations 100 and 380 ppm were applied 
in two separate chambers and maintained throughout this study. 
C0 2 concentrations were monitored and maintained throughout the 
experiments. Plants were grown under a 8/16h light/dark cycle (pho- 
tosynthetic photon flux density 150 umol rar 1 s~') at 21 °C and 70% 
relative humidity. After 4 weeks, the photoperiod was changed to a 
1 6/8 h light/dark cycle for a further 2 weeks. On day 42, samples were 
taken during the middle of the light period and mature expanded 
rosette leaves from 10-15 individual plants were harvested, imme- 
diately frozen in liquid nitrogen, and stored at -80 °C until use. The 
samples were taken from 12 individual pots. 

Morphological data collection 

Scanning electron microscopy and transmission electron micros- 
copy were used to observe the changes of ultrastructure by low C0 2 . 
The number of stomata was counted in four fields of view from the 
fully expanded leaves of no less than eight individual plants for each 
treatment (Supplmentary Fig. SI available at JXB online). 

RNA preparation and sequencing 

Total RNA was prepared with TRIzol (Invitrogen Life Technologies, 
Shanghai, China), according to the manufacturer's instructions. 
Following extraction, total RNA was purified using a RNeasy 
Mini Kit including on-column DNase digestion (Qiagen, Shanghai, 
China). Purified RNA was checked for integrity and quality using 
an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, 
CA, USA). The cDNA library was constructed for sequencing as 
described in Illumina TruSeqTM RNA sample preparation version 
2 guide (catalog no. RS-930-1021). Sequencing was performed using 
a Illumina HiSeq 2000 (Illumina, San Diego, USA). 

Mapping and quantification of sequence reads 

Clean reads were mapped onto the latest A. thaliana Col-0 genome 
assembly (TAIR 10) or a minimal set of coding sequences of the 
TAIR 9 genome release (Gowik et al, 2011) using the bowtie ver- 
sion 0.12.7 (Langmead et al, 2009). The best hit of each read with 
a maximum of three nucleotide mismatches was used (-v 3 —best). 



The raw digital gene expression counts were normalized using 
the RPKM (reads/kb/million) method (Mortazavi et al, 2008; 
Nagalakshmi et al, 2008; Supplementary Tables SI and S2 available 
at JXB online). 

To identify differentially expressed genes, an expression profile 
matrix was built which integrated the digital gene expression count 
for each gene in each library, total gene count for each condition 
were used as background to check if a gene is significantly differen- 
tially expressed in low and C0 2 normal conditions by applying the 
chi-squares test. A FDR-corrected P-value was calculated using 

the formula q(i) = ^-^-C(N) where i represents the ascending 
i 

order of P- values, p(i) represents the fth P-value, C represents a 
chosen constant, and N represents the size of dataset (Benjamini 
and Hochberg, 1995). Significantly differentially expressed 
genes were picked following the criteria P<0.001, FDR<0.025, 
|log 2 Ratio|>1.2. 

Results 

Effects of long-term low C0 2 on biomass growth, 
stomata density, and chloroplast ultrastructure 

C0 2 is the major source of carbon for photosynthesis and 
plays a vital role in plant growth. High C0 2 often increases 
the growth and reproduction of C 3 annuals, whereas low C0 2 
decreases growth (Ward et al, 2000; Ward, 2005). Previous 
studies showed that minimum C0 2 concentrations between 
180 and 200 ppm during the Last Glacial Maximum were 
already stressful on modern C 3 plants (Dippery et al, 1995; 
Ward, 2005); therefore, this work set low C0 2 concentration 
as 100 ppm. Arabidopsis plants grown at 100 ppm for 6 weeks 
were much smaller than those grown under normal C0 2 (380 
ppm) (Fig. 1). In addition, low C0 2 led to a slight delay in 
flowering time (data not shown). The results showed that low 
C0 2 (100 ppm) had a dramatic impact on the growth of the 
C 3 plant Arabidopsis. 

Stomata control the entry of C0 2 into the leaves of plants for 
photosynthesis. There is a strong inverse correlation between 
atmospheric C0 2 and stomatal density (the number of sto- 
mata per unit area) (Franks et al, 2012). This work examined 
the stomatal density of abaxial (lower) leaf blade epidermis 
of Arabidopsis plants grown at either low C0 2 or normal C0 2 
for 6 weeks (Supplmentary Fig. SI). As expected, stomatal 
density was significantly higher (mean±SE 509 ± 59mm -2 ) in 
plants grown at low C0 2 compared to plants at normal C0 2 
(297 ± 54mm- 2 ) (Fig. 2). 

In plants, photosynthesis occurs exclusively in the chlo- 
roplast, and the photosystems (PSI and PSII) exist on the 
thylakoid membrane inside a chloroplast. PSII is limited to 
granal thylakoids, while PSI exists exclusively in the thyla- 
koids exposed to the stroma (Albertsson, 1995; Dekker and 
Boekema, 2005; Sakamoto et al, 2008). The ultrastructure 
of mature leaves under low C0 2 were examined using trans- 
mission electron microscopy, and the size and the arrange- 
ment of bundle sheath cells and mesophyll cells was not 
changed, while Arabidopsis grown under low C0 2 showed 
decreased stacking in chloroplast grana in both mesophyll 
and bundle sheath cells under low C0 2 compared to normal 
CQ 2 (Fig. 3). 
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Fig. 2. Effect of low atmospheric C0 2 on stomatal density. Representative scanning electron micrographs of abaxial (lower) leaf blade epidermis of 
Arabidopsis grown under low C0 2 (1 00 ppm) or normal C0 2 (380 ppm) for 6 weeks. Dashed lines indicate stomata. Bars, 20 |im. 



Some C 4 -cycle genes were upregulated 
under low C0 2 

The mRNA-seq analysis to compare transcriptomes between 
closely related C 4 and C 3 species within the genus Flaveria 
and Cleome using Arabidopsis as the reference genome 
defined a list of enzymes, transporters, and regulatory pro- 
teins required for the core C 4 cycle (Brautigam et al, 2011; 
Gowik et al, 2011). It has been reported that Arabidopsis 
shows the characteristics of C 4 photosynthesis in midveins 
(Brown et al, 2010), but nothing is known about the plastic- 
ity of these characteristics. 



In order to check whether C 4 -related characteristics can 
be regulated by low C0 2 stress, the transcript abundances 
of putative C 4 -related genes were examined. The transcript 
encoding the enzyme phosphoenolpyruvate carboxylase 
(PEPC, At2g42600) showed 2.10-fold higher transcript abun- 
dance, followed by PEPC kinase (PEPC-K, Atlg08650) with 
a 1 .99-fold increase in abundance (Table 1 and Supplementary 
Table S3 available at JXB online). In addition, the transcript 
abundances for the genes encoding alanine aminotransferase 
(Atlgl7290), chloroplast NAD-dependent malate dehy- 
drogenase (At3g47520), pyruvate orthophosphate dikinase 
regulatory protein (At4g21210), inorganic pyrophosphatase 




Table 1 . Transcription abundance of C 4 -cycle genes and C 4 -related transporters 

Reads were mapped onto the latest Arabidopsis thaliana Col-0 genome assembly (gene mapping) or a minimal set of coding sequences of 
the TAIR 9 genome release (core set mapping) using bowtie. Low: low C0 2 , 100 ppm; Nor: normal C0 2 , 380 ppm; rpkm, reads per kilobase 
per million mapped reads. AlaAT, alanine aminotransferase; AspAT, aspartate amino transferase; cpNAD-MDH, chloroplast NAD-dependent 
malate dehydrogenase; Dit, chloroplast dicarboxylate transporter; PEPC, phosphoenolpyruvate carboxylase; PEPC-K, PEPC kinase; 
PEP-CK, PEP carboxykinase; PPA2, inorganic pyrophosphatase 2; PPDK-RP, pyruvate orthophosphate dikinase regulatory protein; PPT1 , 
phosphoenolpyruvate/phosphate translocator 1 ; TPT, triose phosphate transporter; -, no expression detected. 



Gene ID 


Protein 


Gene mapping 






Core set mapping 






Low (rpkm) 


Nor (rpkm) 


Fold change 


Low (rpkm) 


Nor (rpkm) 


Fold change 


At1g08650 


PEPC-K 


26.220 


13.206 


1.985 


30.156 


15.230 


1.980 


At1g17290 


AlaAT 


56.839 


48.904 


1.162 


65.884 


57.025 


1.155 


At1g62800 


AspAT 


0.707 


2.388 


0.296 


0.916 


2.701 


0.339 


At2g 18230 


PPA2 


10.288 


6.468 


1.591 


11.812 


7.459 


1.584 


At2g42600 


PEPC 


197.604 


93.979 


2.103 


227.099 


108.523 


2.093 


At3g47520 


cpNAD-MDH 


74.044 


64.394 


1.150 


85.015 


74.261 


1.145 


At4g21210 


PPDK-RP 


213.075 


191.911 


1.110 


244.822 


221.204 


1.107 


At4g37870 


PEP-CK 


26.606 


43.413 


0.613 


30.589 


50.065 


0.611 


At5g1 2860 


Dftl 


356.959 


275.153 


1.297 


409.991 


317.354 


1.292 


At5g33320 


PPT1 








44.181 


50.558 


0.874 


At5g46110 


TPT 


508.799 


533.917 


0.953 


584.330 


615.779 


0.949 


At5g64280 


Dit2 


27.003 


20.059 


1.346 
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2 (At2gl8230), chloroplast dicarboxylate transporter 1 
(At5gl2860) and 2 (At5g64280) showed trends of upregula- 
tion but their fold changes were less than 2. 

Photorespiratory genes showed trends of upregulation 
under low C0 2 

Low atmospheric C0 2 concentration would increase photores- 
piration, so this work also examined the transcript abundances 
of photorespiration genes. Nearly all genes showed trends of 
upregulation in plants grown under low C0 2 compared with 
those under normal C0 2 (Table 2 and Supplementary Table 
S4 available at JXB online), except for the gene encoding gly- 
cine decarboxylase L-protein (mtLPDl; Atlg48030); however, 
the fold changes were all less than 2. The differential responses 
of genes involved in the photosynthetic light reactions, Calvin 
Benson cycle, and ABA and IAA metabolisms were shown in 
Supplementary Tables S8-1 1 available at JXB online. 

Chloroplast biogenesis- and maintenance-related 
genes showed differential expression in low C0 2 

Given the differential expression of genes involved in chlo- 
roplast biogenesis and maintenance between the C 3 and C 4 
Flaveria species (Gowik etal. , 201 1) and the altered chloroplast 
ultrastructure between low C0 2 and normal C0 2 (Fig. 3), this 
work examined the transcript abundances of genes involved 
in chloroplast biogenesis and maintenance under low C0 2 
and compared them with previously identified genes differ- 
entially expressed between C 3 and C 4 species (Gowik et al, 
2011) (Table 3 and Supplementary Table S5 available at JXB 
online). All the chloroplast biogenesis- and maintenance- 
related genes upregulated by low C0 2 shown in Table 3 were 
also upregulated in C 4 Flaveria species, and five genes down- 
regulated by low C0 2 (At5g52540, Atlg52290, At5g20720, 
At2g32180, and At3gl9820) were also downregulated in C 4 



Flaveria species (Gowik et al, 2011); however, only At44446, 
At5g52540, At3gl7040, and Atlg52290 showed a ratio of 
expression abundance greater than 2. 

Of the genes showing a fold change more than 2, three 
(Atlg44446, At3gl7040, and At5g52540) were enriched in 
C 4 Flaveria species compared to C 3 species. PSII concentra- 
tions are well correlated with chlorophyll b synthesis (Bailey 
et al., 2001), and chlorophyllide a oxygenase (Atlg44446) 
is considered a critical enzyme responsible for chlorophyll 
b synthesis (Yamasato et al, 2005). HCF107 (At3gl7040) 
is a sequence-specific RNA-binding protein and remodels 
local RNA structure in a manner that accounts for its abil- 
ity to enhance translation (Sane et al, 2005; Hammani et al, 
2012). The hcfl07 mutation in Arabidopsis leads to a defec- 
tive PSII (Felder et al, 2001). Although many chloroplast- 
targeted DnaJ proteins have not been characterized, it has 
been hypothesized that chloroplast-targeted DnaJ proteins 
participate in protein folding, unfolding, and assembly pro- 
cesses, and some DnaJ proteins are involved in the stabiliza- 
tion of thylakoid membrane complexes such as photosystem 
II (Chen et al, 2010). Therefore, these three downregulated 
genes were related to reduced PSII and this is in agreement 
with the ultrastructural analysis (Fig. 3). 

Differentially expressed transcription factors 

Ten differentially expressed transcription factors were identi- 
fied (|log 2 Ratio|>1.2) (Table 4). Of these, GOLDEN2-LIKE2 
(GLK2, At5g44190), of the GLK family which is involved in 
chloroplast development (Langdale, 2011), was significantly 
downregulated under low C0 2 . The GLK2 counterpart GLK1 
(At2g20570) was also downregulated in low C0 2 but to a 
lesser extent. 

Waters et al. (2009) identified 20 most upregulated genes 
by GLK1 and GLK2 induction using inducible gene expres- 
sion combined with transcriptome analysis. The current 



Table 2. Transcription abundance of photorespiration genes 

The genes in bold represent these that plays a major function in photorespiration and the knockout results in a low C0 2 -sensitive phenotype 
(Bauwe, 2011). Reads were mapped onto the latest Arabidopsis thaliana Col-0 genome assembly (gene mapping) or a minimal set of coding 
sequences of the TAIR 9 genome release (core set mapping) using bowtie. Low: low C0 2 , 1 00 ppm; Nor: normal C0 2 , 380 ppm; rpkm, reads 
per kilobase per million mapped reads. 



Gene ID 


Enzyme 


Gene 


Gene mapping 




Core set mapping 




Low (rpkm) 


Nor (rpkm) 


Fold change 


Low (rpkm) 


Nor (rpkm) 


Fold change 


At1g11860 


Glycine decarboxylase T-protein 


GLDT1 


909.564 


789.090 


1.153 


1044.524 


909.989 


1.148 


At1g23310 


Glutamate:glyoxylate aminotransferase 


GGT1 


461 .752 


399.672 


1.155 


538.428 


465.798 


1.156 


At1g48030 


Glycine decarboxylase L-protein 


mtLPDl 


250.934 


270.703 


0.927 








At1g68010 


Hydroxypyruvate reductases 


HPR1 


302.977 


286.059 


1.059 


347.873 


329.958 


1.054 


At1g70580 


Glutamate:glyoxylate aminotransferase 


GGT2 


39.407 


23.536 


1.674 








At1g80380 


L-Glycerate 3-kinase 


GLYK 


130.215 


118.051 


1.103 


149.539 


136.138 


1.098 


At2g13360 


Alanine:glyoxylate aminotransferase 


AGT1 


1069.682 


953.983 


1.121 


1228.446 


1100.193 


1.117 


At2g26080 


Glycine decarboxylase P-protein 


GLDP2 


87.272 


75.611 


1.154 


216.751 


173.881 


1.247 


At3g14415 


Glycolate oxidase 


GOX2 


497.994 


473.499 


1.052 


571.960 


546.046 


1.047 


At3g14420 


Glycolate oxidase 


GOX1 


698.380 


617.984 


1.130 


801 .954 


712.742 


1.125 


At4g33010 


Glycine decarboxylase P-protein 


GLDP1 


800.116 


599.088 


1.336 








At4g37930 


Serine hydroxymethyltransferase 


SHM1 


1188.138 


842.873 


1.410 


1364.293 


972.136 


1.403 
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Table 3. Transcript abundance of genes related to chloroplast biogenesis and maintenance 

Reads were mapped onto the latest Arabidopsis thaliana Col-0 genome assembly (gene mapping) or a minimal set of coding sequences of the 
TAIR 9 genome release (core set mapping) using bowtie. Low: low C0 2 , 100 ppm; Nor: normal C0 2 , 380 ppm; rpkm, reads per kilobase per 
million mapped reads. 



Gene ID 


Protein 


Gene mapping 






Core set mapping 




Low (rpkm) 


Nor (rpkm) 


Fold change 


Low (rpkm) 


Nor (rpkm) 


Fold < 


At1 g02560 


CLPP5 (nuclear-encoded CLP 


177.486 


153.738 


1.154 


203.821 


177.293 


1.150 




protease 5), protease subunit 














At1 g06430 


FTSH8 (cell-division 


46.648 


33.676 


1.385 


53.561 


38.836 


1.379 




protease ftsH-8) 














At1 g09340 


CRB (chloroplast RNA binding) 


477.334 


377.221 


1.265 


548.065 


435.050 


1.260 


At1g10350 


Putative DnaJ 


5.668 


8.818 


0.643 


6.507 


10.169 


0.640 




heat-shock protein 














At1 g32080 


Putative membrane protein 


277.058 


238.279 


1.163 


318.113 


274.786 


1.158 


At1 g44446 


Chlorophyllide a oxygenase 


16.090 


45.126 


0.357 


18.475 


52.075 


0.355 


At1 g52290 


Protein kinase-like protein 


6.746 


13.767 


0.490 


- 


- 


- 


At1 g55490 


CPN60B (chaperonin 60 beta); 


294.626 


236.862 


1.244 


346.733 


285.802 


1.213 




RuBisCO large subunit-binding 
















protein subunit beta 














At1 g62750 


SC01 (SNOWY COTYLEDON1); 


369.741 


244.662 


1.511 


443.343 


295.894 


1.498 




elongation factor EF-G 














At1g74730 


Unknown protein 


233.905 


183.326 


1.276 


268.981 


211.675 


1.271 


At2g03390 


uvrB/uvrC motif-containing protein 


43.955 


31 .538 


1.394 


50.468 


36.370 


1.388 


At2g30950 


VAR2 (VARIEGATED 2); 


621.194 


440.354 


1.411 


713.485 


507.890 


1.405 




cell-division protease ftsH-2 














At2g32180 


PTAC18(plastid 


19.442 


27.606 


0.704 


22.323 


31 .836 


0.701 




transcriptionally active 18) 














At2g35490 


Putative plastid-lipid- 


108.463 


88.003 


1.232 


124.626 


101.487 


1.228 




associated protein 3 














At2g46100 


Nuclear transport factor 


63.534 


49.599 


1.281 


72.948 


57.419 


1.270 




2 (NTF2) family protein 














At3g1 7040 


HCF107 (high chlorophyll 


9.314 


21.102 


0.441 


10.694 


24.335 


0.439 




fluorescent 107) 














At3g 19820 


DWF1 (DWARF 1) 


77.847 


87.660 


0.888 


89.382 


101.091 


0.884 


At3g24430 


HCF101 (high chlorophyll 


65.051 


49.942 


1.303 


74.690 


57.594 


1.297 




fluorescence 101) 














At4g24190 


SHD (SHEPHERD)/HEAT 


73.310 


64.708 


1.133 


84.173 


74.623 


1.128 




SHOCK PROTEIN 90-7 














At5g12470 


Unknown protein 


52.805 


35.339 


1.494 


60.629 


40.753 


1.488 


At5g20720 


CPN20 (chaperonin 20) 


247.229 


356.642 


0.693 


283.940 


41 1 .333 


0.690 


At5g42270 


VAR1 (VARIEGATED 1); 


413.651 


316.477 


1.307 


475.001 


364.966 


1.301 




cell-division protease ftsH-5 














At5g52540 


Unknown protein 


16.931 


45.877 


0.369 


19.307 


52.863 


0.365 



change 



work assessed the alteration of these 20 primary targets of 
GLK gene action and found nearly that all of them, except 
COR15a (At2g42540) were downregulated (Table 5 and 
Supplementary Table S6 available at JXB online). COR15a 
was significantly induced under low C0 2 instead, possibly 
because COR 15a is an indirect, secondary target of GLK2 
(Waters et al, 2009). 

Stress-induced mutagenesis pathway was changed 
under low C0 2 

It has been shown that DNA double-strand break-dependent 
stress-induced mutagenesis is important to evolution, through 
producing more mutations under stress in Escherichia coli 



(Cirz et al, 2005; Shee et al, 2011; Al Mamun et al, 2012). 
As a severe stress, can low C0 2 induce more mutagenesis in 
natural populations? This work examined the transcriptional 
changes in genes encoding products related to human DNA 
repair proteins and found that genes involved in damage 
sensing (At5g40450, At2g26980, At4g04720), photoreactiva- 
tion (At3gl5620), homologous recombination (At3g48190), 
nucleotide excision repair (At2g36490, At3g02060, 
At5g04560, Atlg52500, At3g28030, At5g45400), and DNA 
polymerases (At4g32700, Atlg67500) were upregulated by 
low C0 2 (Table 6 and Supplementary Table S7 available at 
JXB online). These results suggest that low C0 2 might induce 
a similar mechanism of DNA double-strand break-depend- 
ent stress-induced mutagenesis to promote evolution. 
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Table 4. Differentially expressed transcription factors using Deseq software 

AP2-EREBP, Apetala 2 ethylene-responsive-element-binding proteins; C2H2, C2H2 zinc finger domain; G2-like, golden2-like; SBR SQUAMOSA 
promoter-binding proteins. P<0.001, FDR<0.025, |log 2 Ratio|>1 .2. 



TF family name 



TF locus ID 



Gene name 



Gene description 



Upregulated under low C0 2 
AP2-EREBP 



C2C2-GATA 



At1 g74930 



At4g26150 



AP2-EREBP 
AP2-EREBP 



At4g34410 
At5g05410 



C2H2 



At5g59820 



Downregulated under low C0 2 

C2C2-CO-like 

SBP 



C2C2-CO-like 
G2-like 



At1g49130 
At2g3381 0 



At4g2731 0 
At5g44190 



MADS 



At5g62165 



ORA47 (Octadecanoid derivative- 
responsive AP2/ERF-domain 
transcription factor 47) 
CGA1 (CYTOKININ-RESPONSIVE 
GATA FACTOR1) 



RRTF1 (redox-responsive 
transcription factor 1) 
DREB2A (dehydration-responsive 
element-binding protein 2A) 



ZAT12 



COLS (CONSTANS-LIKE 8) 
SPL3 (SQUAMOSA PROMOTER 
BINDING PROTEIN-LIKE 3) 



BBX28 

GLK2 (Golden2-like 2) 



AGL42 (AGAMOUS-LIKE 42) 



ORA47 is a regulator of jasmonate 
biosynthesis (Pauwels and Goossens, 2008) 

CGA1 was regulated by light, nitrogen, 
cytokinin, and gibberellic acid, and 
modulated nitrogen assimilation, chloroplast 
development, and starch production 
(Bi era/., 2005; Naito etai, 2007; Mara 
and Irish, 2008; Richteref a/., 2010; 
Hudson et al. , 201 1 ); CGA1 play a key role 
in chloroplast development, growth, and 
divison in Arabidopsis (Chiang etai, 2012) 
RTF1 is involved in redox homeostasis under 
high light stress (Khandelwal etai., 2008) 
DREB2A is involved in dehydration- 
responsive gene expression and 
overexpression of an active form of DREB2A 
results in significant stress tolerance 
to dehydration and significant growth 
retardation (Sakuma et al. , 2006) 
Zat1 2 plays a central role in reactive oxygen 
and abiotic stress signalling in Arabidopsis 
and overexpression of Zat12 in Arabidopsis 
results in the enhanced expression of 
oxidative- and light stress- response 
transcripts (Davletova ef al. , 2005) 

Zinc finger (B-box type) family protein 
SPL3 is involved in regulation of flowering 
and vegetative phase change (Cardon ef al. , 
1997; Wu and Poethig, 2006; Yamaguohi 
ef al. , 2009) 

Zinc finger (B-box type) family protein 
GLK2 is required for normal chloroplast 
development (Fitter ef al. , 2002); GLK2 
together with GLK1 optimize photosynthetic 
capacity by integrating responses to variable 
enironmental and endogenous cues (Waters 
etai., 2009) 

AGL42 is involved in the floral transition and 
RNAi-directed downregulation of AGL24 
results in late flowering (Yu ef al. , 2002) 



Discussion 

This study, as far as is known for the first time, investigated 
responses to low C0 2 at the transcriptome level in model 
plant Arabidopsis. Here, the observed changes of transcrip- 
tomics under low C0 2 are briefly discussed, with particular 
reference to their potential significance for fitness of C 3 plants 
under low C0 2 and potential linkage to C 4 photosynthesis 
evolution. 



Low C0 2 reduced productivity 

Arabidopsis plants grown under low C0 2 had extremely 
small stature compared with plants grown under normal 
C0 2 (Fig. 1). This result is in accordance with previous 
studies on the effect of low C0 2 on plant growth (Ward, 
2005). Arabidopsis grown under low C0 2 has about a 7-day 
delay in flowering time. This has also been observed ear- 
lier (Ward and Strain, 1997) and could be interpreted as 
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Table 5. Transcript abundance of GLK-regulated genes 

Reads were mapped onto the latest Arabidopsis thaliana Col-0 genome assembly (gene mapping) or a minimal set of coding sequences of the 
TAIR 9 genome release (core set mapping) using bowtie. The most upregulated genes by GLK1 and GLK2 induction identified by Waters ef al. 
(2009) were examined and nearly all of them were downregulated by low C0 2 , except COR15a (At2g42540). Low: low C0 2 , 100 ppm; Nor: 
normal C0 2 , 380 ppm; rpkm, reads per kilobase per million mapped reads. CAO, chlorophyllide a oxygenase; CHLH, magnesium chelatase; 
COR15a, COLD-REGULATED 15A; GCN5 related, ornithine N-delta-acetyltransferase; GLK1 , Golden2-like 1 ; GLK2, Golden2-like 2; Lhcb, light 
harvesting complex subunit; MRU1 , mto responding up 1 ; PORB, NADPH:protochlorophyllide oxidoreductase B. 



Gene ID Protein Gene mapping Core set mapping 

Low Nor Fold Low Nor Fold 

(rpkm) (rpkm) change (rpkm) (rpkm) change 



At1g15820 


Lhcb6 


2153.237 


3099.576 


0.695 


2472.919 


3575.038 


0.692 


At1 g44446 


CAO 


16.090 


45.126 


0.357 


18.475 


52.075 


0.355 


At1g76100 


Plastocyanin 


132.858 


261.423 


0.508 








At2g05070 


Lhcb2.2 


431 .847 


1425.637 


0.303 








At2g20570 


GLK1 


32.786 


59.366 


0.552 


37.677 


68.503 


0.550 


At2g34430 


Lhcbl .4 


599.406 


1663.045 


0.360 








At2g35260 


Expressed protein 


74.776 


98.821 


0.757 


85.888 


114.002 


0.753 


At2g39030 


GCN5 related 


0.140 


3.243 


0.043 








At2g42220 


Rhodanese-like 
domain-eontaining protein 


220.272 


249.887 


0.881 


253.769 


288.712 


0.879 


At2g42540 


COR1 5a 


193.204 


31.607 


6.113 


268.502 


47.764 


5.621 


At3g08940 


Lhcb4.2 


326.811 


1274.137 


0.256 








At3g27690 


Lhcb2.4 


136.883 


414.191 


0.330 


157.608 


478.391 


0.329 


At3g56940 


Mg-Proto IX ME cyclase 


428.161 


770.890 


0.555 


491 .646 


889.102 


0.553 


At4g27440 


PORB 


400.475 


873.866 


0.458 








At5g 13630 


CHLH 


456.682 


432.350 


1.056 


524.401 


498.633 


1.052 


At5g35490 


MRU1 


7.921 


28.980 


0.273 


9.095 


33.420 


0.272 


At5g44190 


GLK2 


5.722 


26.922 


0.213 


6.570 


31 .047 


0.212 


At5g54270 


Lhcb3 


2240.644 


3734.380 


0.600 


2573.689 


4307.756 


0.597 



a mechanism to allow for greater accumulation of stored 
reserves that could be allocated to reproduction, resulting 
in increased fitness under low C0 2 (Sage and Coleman, 
2001; Ward, 2005). 

Many studies have shown that atmospheric C0 2 concentra- 
tion negatively regulates stomatal density (Woodward, 1987; 
Beerling et al., 2001; Franks and Beerling, 2009; Doheny- 
Adams et al, 2012; Franks et al, 2012). Paleontological 
research has suggested that the long-term decreases in 
atmospheric throughout the entire evolutionary history 
of vascular plants led to the evolution of high densities 
of small stomata in order to attain the highest g cmax values 
required to counter C0 2 'starvation' (Franks and Beerling, 
2009; Franks et al, 2012). Stomata also exhibit short-term 
adaptive responses to atmospheric C0 2 over much shorter 
timescales. For example, A. thaliana Col-0 grown at high 
C0 2 (720 ppm) had reduced stomata density compared with 
those grown at ambient C0 2 (360 ppm) (Lake et al, 2001). 
In the current work, plants grown under low C0 2 developed 
leaves with higher stomatal density (over 60% increase com- 
pared to normal C0 2 ; Fig. 2), suggesting that the plants 
developed a greater g cmax to counteract the C0 2 limitation 
of photosynthesis. These results suggest that low C0 2 is a 
severe stress to C 3 plants and may greatly reduce C 3 plant 
productivity. 



Responses of genes involved in C 4 photosynthesis and 
photorespiration under low C0 2 

In C 4 plants, C0 2 is initially fixed by the enzyme PEPC into 
a C 4 acid and then transported to the site of Rubisco (Hatch, 
1987). The only photosynthetic gene expression patterns 
common to all independently evolved C 4 lineages are upregu- 
lation of PEPC and downregulation of Rubisco in mesophyll 
cells (Sinha and Kellogg, 1996; Langdale, 2011). Arabidopsis 
has four genes encoding PEPC, and AtPPC2 (At2g42600) is 
the only isoform expressed in leaves. Unlike the other three 
PEPCs, the expression of AtPPC2 is stable and has not 
been reported to be regulated by any stress (Sanchez et al, 
2006; Doubnerova and Ryslava, 2011); however, the cur- 
rent work found that AtPPC2 was upregulated by low C0 2 
(Table 1). The regulators of photosynthetic genes are also 
crucial to maintain C 4 photosynthesis: e.g. plant PEPC activ- 
ity is further regulated through reversible phosphorylation 
by PEPC-K (Nimmo, 2003). Transcripts encoding the C 4 - 
specific regulatory factors PEPC-K and pyruvate orthophos- 
phate dikinase regulatory protein were upregulated as well 
(Table 1). However, changes in other C 4 -related genes were 
less, with fold changes of less than 2. 

When grown in low C0 2 , plants would experience relatively 
high levels of flux through the photorespiratory pathway 



because of the competitive reactions of Rubisco oxygenation. 
In this study, a trend of upregulation of the photorespira- 
tory genes was observed in plants grown under low C0 2 
(Table 2), although most of the genes showed a fold change 
of less than 2. The recent study of transcriptome analysis 
using C 3 , C3-C4 intermediate, and C 4 species of Flaveria 
found that transcript abundances for most genes related 
to photorespiration in the C 3 -C 4 intermediate species 
Flaveria ramosissima were even higher than in the C 3 spe- 
cies Flaveria robusta (Gowik et al, 2011), which is indica- 
tive of the importance of the photorespiratory pathway 
during the evolution of C 4 photosynthesis. The different 
subunits of glycine decarboxylase showed altered expres- 
sion, although the fold changes of these subunits were 
about 0.9-1.3. 

Overall, the data from this study suggest that expression 
of PEPC and PEPC-K is increased under low C0 2 , which 
most likely reflects their potential role for refutation of pho- 
torespired C0 2 under low C0 2 (Sage et al, 2012). For most 
of the other C 4 genes, although trends of upregulation were 
observed, the fold changes were less than 2. Although by 
using expression level changes of all genes under two condi- 
tions as background, this work obtained P-values much less 
than 0.01 for many C 4 -related genes, it is likely that lack of 
biological replicates could had potentially led to an overesti- 
mation of the reliability of statistical tests and caused prob- 
lems in identifying significantly changed genes, especially 
when their fold changes were less than 2. Based on these, this 
work cannot state that low C0 2 induced upregulation of C 4 
genes, except for those genes which showed fold changes over 
2 (e.g. PEPC). 

Readjustment of balance between light absorption and 
C0 2 fixation under low C0 2 

These data on chloroplast ultrastructure and transcript 
abundance of genes involved in chloroplast biogenesis and 
maintenance are consistent with the model for long-term 
photosynthetic regulation by GLK proteins (Waters and 
Langdale, 2009). When light is high and atmospheric C0 2 
is limiting, the rate of C0 2 fixation is insufficient to use 
all of the output of the light-harvesting reactions, result- 
ing in an overly reduced photosynthetic electron transport. 
This triggers a decrease of GLK transcription (GLK1 and 
GLK2; Table 5). Since GLK transcription factors directly 
regulate a large suite of genes involved in light-harvesting 
and thylakoid protein complexes, especially those of PSII 
(Waters et al., 2009), the light-harvesting components in 
the thylakoid membrane LHCB2.2 (At2g05070), LHCB4.2 
(At3g08940), Lhcb3 (At5g54270), Lhcb2.4 (At3g27690), 
and Lhcbl.4 (At2g34430) were downregulated under low 
C0 2 . In addition, the downregulation of the chlorophyl- 
lide a oxygenase gene led to the decrease of chlorophyll b 
synthesis. These results were consistent with the fewer and 
less-stacked grana observed and a higher proportion of non- 
stacked stromal lamellae, as observed in glkl glkl mutants 
(Fig. 3). Therefore, these observed expression changes in 
GLK and the genes regulated by GLK can be interpreted as 
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reflecting the altered balance between C0 2 fixation and light 
absorption. 

Evolutionary implications of plants of to low C0 2 

Growing evidence suggests that all of the basic elements of 
C 4 photosynthesis already existed in C 3 plants. For exam- 
ple, all of the enzymes involved in C 4 photosynthesis exist 
in C 3 plants and play different roles in C 3 plant metabolism 
(Aubry et al, 2011). Some elements controlling the cell spe- 
cific expression of C 4 -related enzymes have been found in C 3 
plants (Brown et al., 2011). Moreover, typical C 3 plants (e.g. 
tobacco and Arabidopsis) show the characteristics of C 4 pho- 
tosynthesis in midveins (Hibberd and Quick, 2002; Brown 
et al, 2010). 

Can some features related to C 4 photosynthesis be 
enhanced under some conditions in a C 3 plant? This work 
showed that under low atmospheric C0 2 , A. thaliana Col-0 
adjusted a series of biological processes, especially the upreg- 
ulation of PEPC and PEPC-K gene expression, and also 
the altered expression of some transcription factors related 
to photosynthesis development, and the downregulation of 
light-harvesting and thylakoid protein complexes. Although 
this study also observed that the majority of the other C 4 - 
cycle genes were upregulated under low C0 2 in Arabidopsis, 
their fold changes were less than 2 and therefore no firm state- 
ments regarding their changes can be made. 

Experiments with more biological replicates and 
Arabidopsis accessions are still needed to firmly conclude 
whether low C0 2 can induce upregulation of other C 4 -related 
genes. Therefore, the results from this paper do not support 
a scenario where low C0 2 acts as a signal to induce C 4 bio- 
chemical features in C 3 plants. It is most likely that the upreg- 
ulation of PEPC and PEPC-K might be a mechanism that 
C 3 plants used to refix photorespired and respired C0 2 and 
also to recapture the released ammonium from photorespira- 
tion and hence increase the competitive advantages under low 
C0 2 conditions. 
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